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Abstract 

This paper introduces a five-parameter lifetime model with increasing, decreasing, up¬ 
side -down bathtub and bathtub shaped failure rate called as the McDonald Gompertz 
(McG) distribution. This new distribution extend the Gompertz, generalized Gompertz, 
generalized exponential, beta Gompertz and Kumaraswamy Gompertz distributions, among 
several other models. We obtain several properties of the McG distribution including mo¬ 
ments, entropies, quantile and generating functions. We provide the density function of 
the order statistics and their moments. The parameter estimation is based on the usual 
maximum likelihood approach. We also provide the observed information matrix and dis¬ 
cuss inferences issues. In the end, the flexibility and usefulness of the new distribution is 
illustrated by means of application to two real data sets. 
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1 Introduction 


The Gompertz (G) distribution, generalizing exponential (E) distribution, is a popular distri- 
bution tha t has b een commonl y used in many applied p roblems for modeling data in biology 


EconomosI (1198211. ge r onto 


ogy 


Brown and Eorbes 




engineering and marketing studies 


Bemmaor and Gladvi (|2012l b A signihcant progress has been made towards the generalization 


and construc tion flexible dis t ributi ons to facilitate better modeling of well-known lifetime data. 
The book by 


Johnson et al. 


(jl99f)l ) provides some applications of the G distribution. 
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In recent years, many authors have proposed distributions which can arise as special sub¬ 
models wit hin the McDona l d gen erated or generalized beta generated (GBG) class of dis¬ 


tributions. 


Alexander et al 


( 20121 1 introduce a class of generalized beta-generated distribu¬ 


tions that have three shape parameters in the generator. They considered eleven differ¬ 
ent parents: normal, log-normal, skewed student-1, Laplace, exponential, Weibull, Gumbel, 
Birnbaum-Saunders, gamma, Pareto and logistic dis tributions. Other generalizations are Mc- 


Donald gamma distribution b y 


Marciano et al 


(j2012l l , McDonald inverted beta d i stribut ion by 


Cordeiro and Lemontel (j2012l l. McDonald nor mal distribution 


Donald extend ed exponential distrib ution by 


distribution by 


Oliveira et al 


(H), 


Cordeiro et al 


McDonald D agum by 




Cordeiro et al 


(I2ni2al l. Mc- 


McDonald half -l ogisti c 
Oluvede and R.aiasoorivJ (2013), 


20121 


McDonald gene ralized beta- b inom ial distribution bv iManoi et al 
3ution by 


istri 


Tahir et al 


3 


(|2013l), McDonald log-logistic 


()2014l ). McDonal d arcsine distri 


( 2 OI 4 I I , and McDonald Weibull distribution by 


3ution by 


Cordeiro et al 


Cordeiro and Lemonte 


((2014) . One of the advantages of 


the McDonald generated distribution lies in its ability of fitting skewed d ata such as other well- 


known distributions in the literature 


(120071 1. 


Mudholkar and NataraianI (|200 


3 ); 


Mudholkar and Wang 


In this paper, we introduce a new five-parameter model called the McDonald Gompertz 
(McG) distribution that includes as special sub-models some recent distributions in the liter¬ 
ature. This distribution offers a more flexible distribution for modeling lifetime data in terms 
of its hazard rate shapes that are decreasing, increasing, upside-down bathtub and bathtub 
shaped. Several mathematical properties of this new model in order to attract wider applica¬ 
tions in reliability, engineering and in other areas of research are provided. 

This paper is organized as follows. In Section [2l we introduce the McG distribution, density 
and hazard functions. Some special models of the new distribution are described in this section. 
In Section[3l we present useful expansions and properties of the cumulative distribution function 
(cdf), probability density function (pdf), kth. moment and moment generating function of the 
McG distribution. Moreover, order statistics and their moments, entropy and quantile measures 
are provided in this section. Estimation of the McG parameters by maximum likelihood (ML) 
method is described in Section (K Finally, application of the McG model using two real data 
sets are considered in Section [5l 




































































2 The McG model 


The generalized b eta distribu t ion o f the first kind (or beta type I) or McDonald distribution 


was introduced by 


McDonaldl (jl984l ). The cdf of the McDonald distribution is given by 


F{x) = I{x^\a/c,b), 0 < X < 1, a, 6 , c>0, 

where I{y; a, b) = Jq — w)^~^dw is the incomplete beta function ratio of 

type I and B{a, b) = — w)^~^dw is the beta function. 

The cdf of McG model can be defined by 


F{y;a,b,c,e,j) = /([I - exp( - -(e^^ - 1))]''; a/c, 6), y > 0, 

7 

where 6,'y > 0. The pdf corresponding to ()2.ip is given by 

cOe^y 9 9 

f{y,a,b,c,9,-f) = exp( -- 1))[1 - exp( -- 1))]“ 

x{l - [1 - exp( - -{e'^y - 1))]"}^-^ 

7 


( 2 . 1 ) 


( 2 . 2 ) 


Here after, we denote a random variable Y with pdf in (12.2p by McG{a, b, c, 9, 7 ). Indeed, the 
McG distribution belongs to McDonald-generalized class of distributions with cdf and pdf as 


1 {y) 

F{y) = I{G‘'{y);a/c,b) = ——— / - w)^~^dw, 

B{a/c,b) Jq 


and 


fiv) = 


c/^.Wb—l 


y(y)G“-^(y)(l-G%)) 


B{ajc, b) ■ 

respectively. The cdf in ()2.ip can be expressed in terms of the hypergeometric function as 

F{y,a,b,c,9,-f) = ^ 2-^1 (a/c, 1 - 6 ; a/c + 1, G^(y)), 

aB[ajc, b) 

where 2 Fi{a, b; c, x) = ^ ascending factorial, and G(y) = 1 — exp( — — 1)) 

is the cdf of G distribution. 


Theorem 2.1. Let f{y,a,b,c,9,^) be the pdf of McG distribution given by (|2.2p . The limiting 
behavior of f{y, a, b,c, 9,'y) for different values of its parameters is given below: 

i. Ifa = 1, then limy^o+ fiy]a,b,c,9,'y) = sirfETj- 

ii. // a > 1, then limj^^Q+ f{y, a, b, c, 9, 7 ) = 0. 

Hi. If a < 1, then limy^Q+ f{y, a, b, c, 9, 7 ) = 00 . 
iv. linij^^+oo f{y] a, b, c, 9, 7 ) = 0. 













Proof. The parts (i)-(iii) are obviously proved. For part (iv), we have 


0 < [1 - [1 - exp{--(e^^ - 1)}]"]^-^ < 1 ^ 

7 

0 < /(y; a, 6, c, 0 , 7 ) < cOe^^ exp(- — 1))[1 — exp(—— — l))]““^/i?(a/c, b). 

7 7 

It can be easily shown that 

lim c0e'^^exp(- — 1))[1 — exp(—— (e'^^ — 1))]““^ = 0, 

y^oo ^ 'y 


and the proof is completed. 
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Density 


Density 


Density 






Figure 1: pdf and hrf of McG model for some values of parameters. 


From equations (j2.1l) and (12.21) . it is easy to verify that the hazard rate function (hrf) of 
the McG distribution is given by 


h{y]a,b,c,e,'y) = 


cOe^'^ exp{—^{e^y — 1 )) 


B{a/c, b) - b) 


9 


x[l - exp(--(e'^^ - 1 ))]“ ^[1 - (1 - exp(--(e^ 2 '- l)))'"]^ ^ y > 0 , 
7 7 


and the corresponding reversed hazard rate function reduces to 


c 6 (e'’'^ exp(—— 1 )) 0 

r{y-a,b,c,9,j) = j, - " [1 - exp(--(e"^ - 1 ))]°-^ 




























































y > 0. 


e 


ci6—1 


x[l-(l-exp(--(e^^-l))r] 

Figured] illustrates some of the possible shapes of density and hazard functions for selected 
values of parameters. For instance, these plots show the hazard function of the new model is 
much more flexible than the beta Gompertz (BG) and G distributions. The hazard rate func¬ 
tion can be bathtub shaped, monotonically increasing or decreasing and upside-down bathtub 
shaped depending on the parameter values. 


BG 


Jafari et al. 

(2014 

), and the beta generalized exponential (BGE) or McE 

Barreto-Souza et al. 


201(1 ) distribu tions for c = 1, a = c and 


ponen tial (BE) iNadaraiah and Kotzl pOOa), the ge neralized Gompertz (GG ) 


7^0, respectively. It also cont ains the beta ex- 


(I 2 OI 3 II. and the Kumaraswam y exponential (KumE) 


GE 


Nadaraiah et al. 


El-Goharv et al 


(I 2 OI 2 II distributions. The 


Gupta and Kundij (119991 ). the G and E distributions are also sub-models. The classes of 


distributions that are included as special sub-models of the McG distribution are displayed in 
Eigure[2| 



Eigure 2: Relationships of the McG sub-models. 


If the random variable Y has the McG distribution, then it has the following properties: 

1. The random variable V = [1 — exp(—— 1))]'^ satisfies the beta distribution with 
parameters a/c and b. Therefore , the ra ndom variable T = — 1) has the BGE (or McG) 


distribution 


Barreto-Souza et al. 


(2010|). Eurthermore, the random variable Y = G = 
























































^ log[l — ^ log(l — follows McG distribution. This result helps us in simulating data from 

McG distribution. The plots comparing the exact McG density function and the histogram 
from a simulated data set with size 100 for some parameter values are given in Figure [3] (left). 
Also, the plots of empirical distribution function and exact distribution function are given in 
Figure [3] (right). These plots indicate that the simulated values are consistent with the McG 
distribution. 

2. If a = z and b = n — i + 1, where i and n are positive integer values, then the F{y; a, b, c, 6, 7 ) 
is the cdf of the zth order statistic of GG distribution. 

Histogram Empirical Distribution 




Figure 3: The histogram of a generated data set with size 100 and the exact McG density (left) 
and the empirical distribution function and exact distribution function (right). 


3 General properties 

In this section, some properties of McG distribution are considered. 


3.1 A useful expansion 


We derive some expansions for the cdf, kth moment and moment generating function of the 
McG distribution. The binomial series expansion is defined by 


{1-zr 


00 





00 


J=0 


F(m + 1) 
T{m-j + 1) j\ ’ 


(3.1) 


where |z| < 1 and m is a positive real non-integer. 

The following proposition reveals that the McG distribution can be expressed as a mixture 
of distribution function of GG distribution, whereas Proposition [22] provides a useful expansion 





















for the pdf in ()2.2I1 . 

Proposition 3.1. The cdf in m is a mixture of distribution function of GG distribution on 
the form 

oo oo 

F{y;a,b,c,9,'y) = = '^pjGj{y), (3.2) 

j=0 j=0 

where pj = B(a/c,b)r{L^j)’jl{a/c+j) is the distribution function of a 

random variable which has a GG distribution with parameters 9, 7 and a + jc. 

using binomial expansion (13.ip to the term [G{y)]°‘^^^ in (j3.2l) . we have 

10(9)1“+'“ = (1-0(9))'= 

fc =0 ^ 

= fE(-l)'=(“V')(')[0(9)|0 

Now, (j3.2h becomes 

OOCXDCXD /|•\/7\ 

F{y;a,b,c,9,'y) = E E Eft(-i)‘'+' (“ fj (J [g(9)1’- = E ^[g(9)1“, 

j=0 r=0 k=r ^ / \ / j,_q 

where 6, = Ef-o Sr.rf>r(-l)'=+’'(“t'“) (S' 

Proposition 3.2. The pdf of McG can be expressed as an infinite mixture of GG densities 
with parameters 9, 7 and {a + jc) given by 

00 00 

f{y,a,b,c,9,j) = jc)c/(?/)[G(y)]“+^''“^ = '^PjOjiv), 

j=0 r=0 

where gj{y) = (a + jc)g{y)[G{y)]°‘~^^'^~^. We can write the pdf of McG as 

00 

f{y,a,b,c,9,j) = g{y) Eo10(9)I“, 

r =0 

where e. = E”. (-l)'+*+ 7,Vca«,-r) GW) C) ■ 

3.2 Moments and generating function 

In this section, we deal with the basic statistical properties of McG distribution such as the 
A)-th moment and generating function in the following propositions. 

Proposition 3.3. The k-th moment of McG distribution can he expressed as a infinite mixture 
of the k-th moment of GG distributions as follows: 

00 00 

j=0 j=0 







where 


E^r) =E E (“ ^ r') +'>!’ 


-1 


ifc+l 


7 (r + 1)' 


i=0 r=0 

and Ujk = 0{a + jc)T{k + 1). 

Proposition 3.4. An explicit expression for the moment generating function of McG distri¬ 
bution follows from Provosition \3.S\. 


where 


poo 

Mrit) = e^^J2Pj{a + jc)g{y)[G{y)r+^^-^ = (t), 


MvAt) = 


j=0 
(o + jc)9 


j=0 


7 


EE(-i)' 


i=0 k=0 


a + jc—l\ft/^\ r(/c + l) 


k ) f(^±Z£)®]fc+i ■ 


3.3 Order statistics 


Order statistics make their appearance in many areas of statistical theory and practice. Let 
the random variable Yi-n be the iih. order statistic {Yi-n < Y 2 :n < • • • < Yn-.n) in a sample of 
size n from the McG distribution. The pdf and cdf of Yi-n for z = 1,2,..., n are given by 

1 


fi:n{y) — 


B{i, n — i + 1) 
1 


f{y)[F{y)r^[l-F{y)] 




S(i,n-i + l)^V k 


n — i 


(-i)7(!/)[n!/)i 


fc+2 —1 


and 


F 


■.n{y) = f fi-.n{t)dt 

Jo 


B{i,n — i + 1) k + i \ k 




(3.3) 


(3.4) 


respec tively, where F{y) = FG{y). We use throughout an equation by (iGradshtevn and R.vzhik 

(120071 ). page 17) for a power series raised to a positive integer m given by 


I ='^Cm,rU^, (3.5) 

V ?’=0 / r =0 

where the coefficients Cm,r (for r = 1, 2,...) are easily determined from the recurrence equation 

r 

c-m,r — {j" ho) ^ ^ \k (m T 1) r T fc] Cm,i —ki 

k=l 

where Cm,o = b^f. Hence, the coefficients Cm,r can be calculated from Cm, 0 ) • • • j Cm,r-i and 
therefore, from the quantities bo,---,br. Using (j3.5|) . the equations ([33]) and (1331) can be 
written as 

^ n—i oo / -x 

1 r , .su(n — i^ 


fi:n{y) — 


B{i, n — i + 1) 


EEr^M) 


fc=0 r=l 


k + F k )^^^^My)[Giy)r\ 















Fr.n{y) = 


EErT7(-‘)‘ 


B{i,n — i + 1) ^ ^ /c + z' 

^ ’ ’ A:=0 r=0 


n — I 


Ci+k,r[G{y)\ 


An explicit expression for the sth moments of Yi-n can be obtained as 


= 


1 


B{i, n — i + 1) 

erjs + i) 

B{i, n — i + \) 


EErT7<-i) 

k=0 r=l 
n—i oo 

EErlrTl-i) 


/c=0 r=l 


n — I 


n — i 


r+oo 


Ci+k,r / t^C/(t)[G(t)]'’ 




21=0 22=0 

3.4 Quantile measures 


^ ^ /rA - 1\ (-l)'^+'2 i(n+i). 6)(ii + 1) ^ -1 

.^n.^nV yr(z2 + l) ^ 7 ^ + ^ ^ 


In this section, we consider the effect of each shape parameters a, b and c on the skewness and 
kurtosis of the McG distribution. To illustrate this effect, we use measures based on quantiles. 
The quantile function of the McG{a,b, c,'y,6) distribution say Q{t) can be obtained as 


Q{t) = ^ log(l - I log(l - 0 < t < 1, 


where Qa/c,b{t) = denote s the fth quantile of b eta 


istribution with parameters 


a/c and b. The Bowley skewness fsee iKennev and Keeping] (jl962l )l based on quantiles can be 
calculated by 


B = 


and the Moors kurtosis (see iMoorsI (jl988l B is defined as 


Q(f)-2 g(i) + g(i) 


,, Q(|)-Q(i) + Q(i)-Q(i) 

Q(|)-Q(i) 

where Q{.) denotes the quantile function. These measures are less sensitive to outliers and 
they exist even for distributions without moments. For the standard normal and the classical 
standard t distributions with 10 degrees of freedom, the Bowley measure is 0. The Moors 
measure for these distributions is 1.2331 and 1.27705, respectively. 

In FigurelH we plot Bowley measure (holding b = 0.5, 7=1 and 9 = 0.1 fixed) as a function 
of c for fixed values of a (left) and Bowley measure (holding o = 0.5, 7 = 1 and 0 = 0.1 fixed) 
as a function of c for some values of h (right). In Figure El we plot Moors measure (for h = 0.5, 
7 = 1 and 0 = 0.1) as a function of c for selected values of a (left) and Moors measure (for 
a = 0.5, 7 = 1 and 0 = 0.1) as a function of c for some values of b (right). These plots indicate 
that these measures can be sensitive to the three shape parameters a, b and c. 
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Figure 4: The Bowleys skewness of the McG distribution as a function of c. 


Moors kurtosis 


Moors kurtosis 




Figure 5: The Moors skewness of the McG distribution as a function of c. 

3.5 Entropy 

The entropy of random variable is defined in terms of its probability distribution and can 


be shown to be a good measure of randomness or uncert ainty. T 


continuous random variable V with pdf f(y) is defined bv IShannonI (jl948l f as 


l e Sha nnon’s entropy of a 


1*00 

J^Shif) = -Ef [log f(Y)] = - f (y) log f{y)dy. 

Jo 

Hence, the Shannon entropy for McG distribution can be expressed in the form 
Eshif) = - ^/l - iE{Y) + 6 »/ 7 My ( 7 ) + (a - l)C(a, b) + {b - 1)C(6, a), (3.7) 


where C(r, s) = ^p{r + s) — ^lJ{r) and V’(-) represents the digamma fun ction. The last two terms in 
(13.71) follows immediately from the first two conditions in Lemma 1 of lZoerafos and Balakrishnan 





















































([20091). The Renyi entropy is defined by 


P 7—00 

where p > 0 and p ^ 1. The Shannon entropy is derived from \mip^iHp[f). An explicit 
expression of Renyi entropy for McG distribution is obtained as 

HpU) = - log(6>) + —^ log( ^ ) + -7— log(R(ap - p + cj + 1, p)) 

1 — p B[a/c,b) 1 —p 

+rb ‘<>e(E(-i)' (*'’ 7)®{(1 - -iio I°e(i - c/))'-'}. 

P j=0 \ J / 

where U has a beta distribution with parameters ap — p + cj — I and p. 


4 Estimation 


Let Yi,... ,Yn he a random sample of size n from the McG{a, b, c, 9, 7 ) distribution and © = 
{a,b,c,9,^) be the unknown parameter vector. The log-likelihood function is given by 

n n ^ 

l{&) = nlog{c9)-n\og{B{a/c,b))+j^yi -- 1) 

i=l i=l 

n n 

+ (a - 1) X] + (& - 1) X] - (1 - U)''), (4.1) 

i=l i=l 

where ti = ex.p{—^{e'^^* ~ !))• The maximum likelihood estimation (MLE) of 0 is obtained by 
solving the nonlinear equations, 17(0) = {Ua{®), Ub{&), Uc{&), Ug{®), U.y{@))'^ = 0, where 


17a (0) 

Ubi®) 

Uci®) 


= n/c[ip{a/c + 6) - V’(a/c)] + ^ log(l - U), 

i=l 
n 

/c + b)- ijjib)] + log(l - (1 - tf)), 
= n/c — nal(?['Ll){alc + b) — V’(a/c)] 


da 

dl{®) 

db 

dl{®) 

dc 


= n\w(a 


.( 5 _ - ti) 


Uei®) = 


dl{®) 

do 


fcr ’ 

n n 

= n/e - 1/7 - 1) + (a - l)/7 Y 


2=1 


2=1 


- 1 ) 

l-ti 


-c{b-l)hY- - 1 


- 1 ) 


2 = 1 




UY®) = 


dl{®) 

d'j 


Yy^~^ ~ -1) 


2=1 


2=1 


















i=l 

n 


+6»(6 - l)c/7^ ^ 


2=1 


l-ti 

ti{l - - 1) 

1 - (1 - uY 


We need the observed information matrix for interval estimation and hypotheses tests on 
the model parameters. The 5x5 Fisher information matrix, J = J„(0), is given by 



Jaa 

Jab 

Jac 

Jad 

Ja^ 


Jba 

Jbb 

Jbc 

Jbe 

Jb'j 

J = - 

Jca 

Jcb 

Jcc 

Jed 

Jcy 


Joa 

Jeb 

Jqc 

Jee 

JO'y 


_ J'ya 

J^b 

J'yc 

J'yd 

J 


where the expressions for the elements of J are 
Jaa = + h)- Y'{a/c)], 

Jac = =-^[Y'ia/c + b)-Y'ia/c)], 


JaO — 


Ja') — 


Jbb — 

Jbc 


dadc 

d^ije) _ 1 A tije^y^ - 1) 
dadd 7 ^ 1 — t j ’ 

52/(0) _ 0 A UijVie'ry^ - + 1) 


— V 

5a57 7^ ^ 

52/(0) 


56^ 

52/(0) na 

565 c c 


l-ti 

= n[Y'{a/c + 6) - Y'{b)], 


^Y'ia/c + b) - 


(1 - ti)''log(l - ti) 


2=1 


1- (1-ti)^ 


^ _ 52/(0) _ c ti{i - uy-Ye^y* - 1) 

^be - o, o/. - 2^ ^ 


Jby — 


5650 7 ^ 1 — (1 — tj)^ ’ 

52 /( 0 ) c0 ^ ti(l - ti)‘^“H72/ie^''^ - e^y^ + 1) 


5657 


ca til, 
'Z2 

’ i=l 


1 - (1 - ti 


52 /( 0 ) n 2na. , , , ,, , / , ^ 

Jcc = o 9 =—-^[V'(a/c + 6) - ^'(a/c)] 


5c2 


c2 C^ 


Jc6 — 


na?. (1 — ti)‘^(log(l — ti ))2 

+—[V' («/c + &) - V’ (a/c)] - (6 - 1) ^- [1 _ (1 _^.)c ]2 ’ 

52/(0) _ 6 - 1 ■A tj(l - tiY~^{e^y^ - l)[clog(l - ti) + 1 - (1 - tiY] 


dcdO 


7 


E- 

2=1 


[1 - (1 - UY? 


^ _ 52/(0) _ 0(6 - 1) ti(7yie^2'* - _t_ i)[clog(l - fi) + 1 - (1 -/i)‘=] 

~ ~ ~ ^,2 A 


dcd'f 


T 


2=1 


(i-t.)i-1i-(i-AF 



5 ^/( 0 ) n a — 1 ^ — 1)^ c(6 — 1) 

^00 o/i9 ~ 7i9 9 / V ' 


502 


02 ^2 ^ 


7 


ii(l - UY 2(e7?/i _ l)2{ctj + (1 - tif - 1} 


(1 - (1 - tiY? 


JQ'y — 


X:(e™ - ,y.e^y. - D - E ^ ^ 


2 = 1 

52/(0) _ j_ 

5057 ^ (1 - - (1 - 

X [(1 - tiYiee^y^ +ta--f-e) + couie'^y^ -1) + 7(1 - u) + 0(1 - 

(o — 1) ti{e^y* — ^Vie^y^ — l){ 0 e^y^ + 7tj — 7 — 0) 


+ 


--yO 


E 

2=1 


Y-k) 


72 


Jii = ^^ = ^E(70*e^"^-e^"*-7Ve^"72 + l) 

^ ^ i=l 


c(6 - 1)02 " YY - tif-^Yvie^y^ - e^y^ + 1)2(2c(1 - tif + (1 - tif - c - 1) 


T 


E 

2=1 

n 


{i-{i-kYY 


^ 2Yyie*y* - 2je'fy^ + 27 + 

7 ^ 1 - (1 - ^5 

+0(e^2/i _ 1)2 _ 26-fyie^yYe'^y^ - 1)] 

_ ^ ^_^3y2^'yy, ^ 2Yyie'^y^ - 27 (e^^» - 1) + O^yie'^y^ - e^y^ + 1)^] 

1 ti 


2=1 


(a - 1)02 ^ tji-fyie'ryi - e^y^ + 1)2 


7 


2=1 


(l-t.)^ 


Under conditions that are fulfilled for parameters in the interior of the parameter space but 
not on the boundary, asymptotically 

V^(0-0)~ Ar5(O,I(0)-'), 


where /(0) is the expected information matrix. This asymptotic behavio r is valid if /(0) 


replac ed by Jn(0) , i-e., the observed information matrix evaluated at © 

umd)- 


Cox and Hinklev 


5 Application of McG to two real data sets 

In this section, two real data sets are considered to illustrate that the McG model can be a good 
lifetime distribution comparing with main three submodels; BG, KumG and McE distributions. 
In both examples, we obtain the MLE and their corresponding standard errors (in parentheses) 
of the model parameters. The model selection is carried out using minus of log-likelihood 
function (—log(L)), Kolmogorov-Smirnov (K-S) statistic with its p-value, Akaike information 























criterion (AIC), Akaike information criterion corrected (AICC), Bayesian information criterion 
(BIC) and likelihood ratio test (LRT) with its p-value. Furthermore, we plot the histogram for 
each data set and the estimated pdf of the four models. Moreover, the plots of empirical cdf 
of the data sets and estimated cdf of four models are displayed. 


Example 5.1. The data set have been obtained from 


of 50 devices. Also, it is analyzed by 


El-Goharv et al. 


Aarsei ll981i) and represent s the lif etimes 


t '2013i) and 


Jafari et al 


(MU)- The 


results whieh are given in Table 1 indicate that the MeE model is not suitable for this data set 
based on K-S statistie. The MeG model has the lowest — log{L), AIG and AIGG values among 
all fitted models, but the BG model has the lowest BIG value among all other fitted models. A 
comparison of the proposed distribution with some of its submodels using p-value of LRT shows 
that the McG model yields a better fit than the other three distributions to this real data set. It 
is also clear from Figure @ that the McG distribution provides a better fit and therefore be one 
of the best models for this data set. 


Example 5.2. The data set represents the strengths of 1.5 em glass fibers, measured at 
the National Physical Laboratory, England. Unfortunately, t h e uni ts of measurement are 


not given in the paper. I t is obtained from 


Barreto-Souza et al. 


Smith and Naulon u98il) and also analyzed by 


(201m)- The K-S statistie indieates that all fitted models are good eandi- 


dates for this data set, but based on all other criteria: — log(L), AIG, AIGG, BIG and p-value 
of LRT, we can infer that the best model is the McG distribution. Similar results are concluded 
from Figure 
























Table 1: MLEs of the model parameters for the data of lifetimes of 50 devices, the corresponding 


SEs and the K-S, AIC, AICC, BIC and LRT statistics. 



Distribution 


BG 

KumG 

McE 

McG 

a 



0.8643 

0.2619 

(s.e.) 

(0.0392) 

(0.0871) 

(0.4621) 

(0.0656) 

b 

0.2467 

0.6063 

0.0706 

0.0752 

(s.e.) 

(0.0448) 

(0.3780) 

(0.0927) 

(0.1029) 

c 

— 


2.7127 

3.7652 

(s.e.) 

— 

(0.0871) 

(3.7438) 

(0.9946) 

e 



0.2826 

0.0012 

(s.e.) 

(0.0001) 

(0.0001) 

(0.3727) 

(0.0001) 

7 


0.0766 

— 

0.0875 

(s.e.) 

WOBIlBliM 

(0.0282) 

— 

(0.0001) 

- log(L) 

220.6714 

221.9666 

237.8158 

219.0041 

K-S 

0.1322 

0.1367 

0.1916 

0.1216 

p-value (K-S) 

0.3456 

0.3072 

0.0507 

0.4509 

AIC 

449.3437 

451.9331 

483.6316 

448.0081 

AICC 

450.2326 

452.8220 

484.5205 

449.3718 

BIC 

456.9918 

459.5812 

491.2797 

457.5682 

LRT 

3.3356 

5.9249 

37.6233 

— 

p-value (LRT) 

0.06779 

0.0149 

0.0000 

— 


Table 2: MLEs of the model parameters for the strengths of 1.5 cm glass fibers data, the 
corresponding SEs and the K-S, AIC, AICC, BIC and LRT statistics. 



Distribution 


BG 

KumG 

McE 

McG 

a 

1.6907 

1.8946 

9.3276 

0.7940 

(s.e.) 

(0.8664) 

(1.2350) 

(2.6891) 

(0.2355) 

b 

27.7434 

4.2814 

93.4655 

0.1248 

(s.e.) 

(87.8726) 

(15.7612) 

(109.2042) 

(0.1786) 

c 

— 

1.8946 

22.6124 

192.1704 

(s.e.) 

— 

(1.2350) 

(14.0674) 

(307.3698) 

e 



0.9227 

0.0009 

(s.e.) 


(0.0509) 

(0.3443) 

(0.0001) 

7 

2.7156 

2.3019 

— 

5.2013 

(s.e.) 

(0.8448) 

(1.6517) 

— 

(0.4018) 

- log(L) 

14.2158 

14.0305 

15.5995 

11.4208 

K-S 

0.1324 

0.1313 

0.1466 

0.1159 

p-value (K-S) 

0.2186 

0.2271 

0.1335 

0.3651 

AIC 

36.4317 

36.0610 

37.2569 

32.8417 

AICC 

37.1213 

36.7506 

37.9466 

33.8943 

BIC 

45.0042 

44.6335 

45.8295 

43.5573 

LRT 

5.5900 

5.2193 

8.3572 

— 

p-value (LRT) 

0.0181 

0.0223 

0.0038 

— 
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Figure 6: Plots of the estimated pdfs and cdfs BG, KumG, McE and McG models using data 
of lifetimes of 50 devices. 

Histogram Empirical Distribution 




X 


X 


Figure 7: Plots of the estimated pdfs and cdfs BG, KumG, McE and McG models using the 
strengths of 1.5 cm glass fibers data. 
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